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Abstract 

We demonstrate that a tensor product structure and optical analogy of quantum entanglement 
can be obtained by introducing pseudorandom phase sequences into classical fields with two orthog¬ 
onal modes. Using the classical analogy, we discuss efficient simulation of several typical quantum 
states, including product state, Bell states, GHZ state, and W state. By performing quadrature 
demodulation scheme, we propose a sequence permutation mechanism to simulate certain quantum 
states and a generalized gate array model to simulate quantum algorithm, such as Shor’s algorithm 
and Grover’s algorithm. The research on classical simulation of quantum states is important, for 
it not only enables potential beyond quantum computation, but also provides useful insights into 
fundamental concepts of quantum mechanics. 
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Introduction 


It has been widely known that fast algorithm is not available for many famous computing 
problems according to classical computational complexity theory [TJ. Quantum computa¬ 
tion, however, enormously promotes computational efficiency with accumulation of compu¬ 
tational complexity no faster than polynomial rate under linear increasement of input size 
by using several basic and purely physical features of quantum mechanics, such as coherent 
superposition, parallelism, entanglement, measurement collapse etc. |2]|. The acceleratory 
ability of quantum computation is related to tensor product and quantum entanglement, of 
which the latter one is essential to realize quantum computation and quantum communica¬ 
tion m- Yet quantum algorithm is difficult to be realized for restrict of quantum system 
controllability, decoherence property and measurement randomness [SHE]. 

Recently, several researches have proposed a new concept of realization of classical en¬ 
tanglement based on classical optical fields by introducing a new freedom degree to realize 
tensor product in quantum entanglement [9]. Further, Ref. ra proposed that phase mod¬ 
ulation by orthogonal pseudo-random sequence is able to simulate quantum entanglement 
effectively. In this scheme, these classical fields with an increased freedom degree not only 
realized tensor product structure but also simulated the nonlocal property of quantum entan¬ 
glement by using property of orthogonal pseudo-random sequence like orthogonality, balance 
and closure. It is interesting that this simulation creates a Hilbert space of n2 n dimensions 
which is larger than quantum mechanics. Different from quantum measurement involving 
collapse, these classical fields are measured in intensity by optical detectors directly. 

The classical simulation of quantum systems, especially of quantum entanglement has 
been under investigation for a long time HM3]. In addition to potential practical applica¬ 
tions in quantum computation, research on classical simulation systems can help understand 
some fundamental concepts in quantum mechanics. However, it has been pointed out by 
several authors that classical simulation of quantum systems exhibits exponentially scaling 
of physical resources with the number of quantum particles [[3|IT3]. In Ref. p3], an optical 
analogy of quantum systems is introduced, in which the number of light beams and opti¬ 
cal components required grows exponentially with the number of qubits. In Ref. HI, a 
classical protocol to efficiently simulate any pure-state quantum computation is presented, 
yet the amount of entanglement involved is restricted. In Ref. jB], it is elucidated that in 
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classical theory, the state space of a composite system is the Cartesian product of subsys¬ 
tems, whereas in quantum theory it is the tensor product. This essential distinction between 
Cartesian and tensor products is precisely the phenomenon of quantum entanglement, and 
viewed as the origin of the limitation of classical simulation of quantum systems. Therefore 
it is of great significance in the classical simulations to realize tensor product PHD]. 

In wireless and optical communications, orthogonal pseudorandom sequences have been 
widely applied to Code Division Multiple Access (CDMA) communication technology as a 
way to distinguish different users [IS HSJ. A set of pseudorandom sequences is generated 
from a shift register guided by a Galois held GF(p), which satisfies orthogonal, closure and 
balance properties [I6j- In Phase Shift Keying (PSK) communication systems [IT], pseu¬ 
dorandom sequences are used to modulate the phase of the electromagnetic/optical wave, 
where a pseudorandom sequence is mapped to a pseudorandom phase sequence (PPS) values 
in (0, 27 r/p, ■ ■ ■, 27t (p — 1) /p}. Guaranteed by the orthogonal property of the PPS, different 
electromagnetic/optical waves can transmit in one communication channel simultaneously 
without crosstalk, and can be easily distinguished by implementing a quadrature demodu¬ 
lation measurement ra- 

In this paper, we propose a new mechanism based on circular demodulation to realize 
simulation of certain quantum state represented by these classical optical fields. Besides, 
classical simulation of some other typical quantum states is discussed, including product 
state, Bell states, GHZ and W states. Furthermore, we propose a new scheme to simulate 
quantum computing by constructing certain quantum states based on an array of several 
mode control gates. Finally, we use this method to simulate Shor’s algorithm |18; [19j and 
Grover’s algorithm m- 

The paper is organized as follows: In Section [T] we introduce some preparing knowledge 
needed later in this paper. In Section |TT] the existence of the tensor product structure in our 
simulation is demonstrated and the classical simulation of several typical quantum states is 


analyzed. In Section III a generalized gate array model to simulate quantum computation 


is proposed. Finally, we summarize our conclusions in Section IV 


3 




I. PREPARING KNOWLEDGE 


In this section, we introduce some notation and basic results needed later in this paper. 
We first introduce pseudorandom sequences and their properties. Then we discuss the 
similarities between classical optical fields and single-particle quantum states. Finally, we 
introduce the scheme of modulation and demodulation on classical optical fields with PPSs. 

A. Pseudorandom sequences and their properties 

As far as we know, orthogonal pseudorandom sequences have been widely applied to 
CDMA communication technology as a way to distinguish different users mm- A set of 
pseudorandom sequences is generated from a shift register guided by a Galois Field GF (p), 
which satisfies orthogonal, closure and balance properties. The orthogonal property ensures 
that sequences of the set are independent and distinguished each other with an excellent 
correlation property. The closure property ensures that any linear combination of the se¬ 
quences remains in the same set. The balance property ensures that the occurrence rate 
of all non-zero-element in every sequence is equal, and the the number of zero-elements is 
exactly one less than the other elements. 

One famous generator of pseudorandom sequences is Linear Feedback Shift Register 
(LFSR), which can produce a maximal period sequence, called m-sequence [IS]. We consider 
an m-sequence of period N — 1 (N = p s ) generated by a primitive polynomial of degree s 
over GF(p). Since the correlation between different shifts of an m-sequence is almost zero, 
they can be used as different codes with their excellent correlation property. In this regard, 
the set of N — 1 m-sequences of length N — 1 can be obtained by cyclic shifting of a single 
m-sequence. 

In this paper, we employ pseudorandom phase sequences (PPSs) with 2-ary phase shift 
modulation, which is a well-known modulation format in wireless and optical communi¬ 
cations, including PSK [IT]. Different from Ref. [10], we choose GF(2) instead of GF(4) 
because the correlation measurements for Bell inequality are not discussed in this paper. 
We first propose a scheme to generate a PPS set 5 = |A®, A^ 1 ),..., over GF(2) 

[IBJ. A^ is an all-0 sequence and other sequences can be generated by using the method as 
follows: 
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(1) given a primitive polynomial of degree s over GF(2), a base sequence of a length 2 s — 1 
is generated by using LFSR; 

(2) other sequences are obtained by cyclic shifting of the base sequence; 

(3) by adding a zero-element to the end of each sequence, the occurrence rates of all 
elements in all sequences are equal with each other; 

(4) mapping the elements of the sequences to {0,7 t/ 2}: 0 mapping 0, 1 mapping tt/2. 
Further, we define a map / : A —> e lX on the set of S, and obtain a new sequence 


set fl = 1 tp l '- 3) 


pU) = e ixU) , j = 0 


In 


Fig- [lj 


we demonstrate the relationship 


between time slots, an m-sequence and phase sequence A*' 4 -* with N phase units A^ = 
[AS 0 , A?, • • • A^ ]. For better understanding our scheme, the PPSs in the cases of modulating 
7 classical optical fields is illustrated below. Using the method mentioned in section |T} an 
m-sequence of length 2 3 — 1 is generated by a primitive polynomial of the lowest degree over 
GF( 2), which is 1110010- Then we obtain a group that includes 8 PPSs of length 
8: {AW,..., AT)|, of which in exception to A^°\ all PPSs can be used to modulate classical 
optical fields to simulate quantum states of up to 7 particles expect A^°y for example, 
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According to the properties of m-sequence, we can obtain following properties of the set 
fl, (1) the closure property: the product of any sequences remains in the same set; (2) the 
balance property: in exception to (p(°\ any sequences of the set satisfy 
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je,Jj) _ 
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Xy( A * )+0 ) = o,V0 e M; 


( 2 ) 


(3) the orthogonal property: any two sequences satisfy the following normalized correlation 
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FIG. 1: The time sequence relationship of the PPS is shown. 

1 , i = i 

0 , i^j 

In conclusion, the map / corresponds to the modulation of PPSs of on classical optical 
fields. According to the properties above, the classical optical fields modulated with different 
PPSs become independent and distinguishable. 

B. Similarities between classical optical field and single-particle quantum states 

We note the similarities between Maxwell equation and Schrodinger equation. In fact, 
some properties of quantum information are wave properties, where the wave is not necessary 
to be a quantum wave |J3j. Analogous to quantum states, classical optical fields also obey 
a superposition principle, and can be transformed to any superposition state by unitary 
transformations. Those analogous properties made the simulation of quantum states using 
polarization or transverse modes of classical optical fields possible [2T1 1231 . 

We first consider two orthogonal modes (polarization or transverse), as the classical sim¬ 
ulation of quantum bits (qubits) |0) and |1) [T31I2T] . 

(4) 

Thus, any quantum state of a single particle can be simulated by a corresponding classical 
mode superposition field, as follows 

W) = « |0) + p I 1 ), M 2 + \P\ 2 = 1, {a,P e C). (5) 
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Obviously, all the mode superposition fields can span a Hilbert space, where we can perform 
unitary transformations to transform the mode state. For example, the unitary transforma¬ 
tion U (y, Q) is defined 


U(x,0) 


fj-xWx cos 9+cr y sin 9) 


cosy 


-ie lG sin x 


le sin x cos X 


( 6 ) 


where a x , a y are Pauli matrices. The modes |0) and |1) can be transformed to mode 
superposition by using U (y, 9), respectively, as follows 


U(x,0)\0) = cos x |0) + ie l ° sin y 11), (7) 

U( y, 6») |1) = cos y 11) — ie~ ld sin y |0). 


Now, we consider some devices with one input and two outputs, such as beam or mode 


splitters, which split one input field \'ipi n ) — a |0) + /3 |1) into two output fields 
/■ For the case of beam splitters, the output fields are 




(a) 

out 


and 


and 


i’out) = C a (a |0) +/3e^“ |1)) 
^out) = Cb (oc |0) + |1)) with an arbitrary power ratio | C a \ 2 : \Cb\ 2 between the 


output beams, where 0 a ^ are the additional phases due to the splitter. For the case of mode 


splitters, the output fields are 


= OLe l ^ a |0) and ip ^ t ^ = P el<t>b 11), where are also 

the additional phases. Conversely, the devices can act as beam or mode combiners in which 
beams or modes from two inputs are combined into one output. 


C. Modulation and demodulation on classical optical fields with pseudorandom 
phase sequences 

We first consider the modulation process on a classical optical field with a PPS. Similar 
to PSK system, choosing a PPS in the set of S, the phase of the field can be modulated 
by a phase modulator (PM) that controlled by a pseudorandom number generator (PNG), 
the scheme is shown in Fig. [2} If the input is a single-mode field, it can be transformed to 
mode superposition by performing a unitary transformation after the modulation. 

Then we consider the quadrature demodulation process of a modulated classical optical 
field. Quadrature demodulation is a coherent detection process that allows the simultaneous 
measurement of conjugate quadrature components via homodyning the emerging beams with 
the input and reference fields by using a balanced beam splitting, where the reference field is 


7 








PNG 


Input field rr~i— —1 Output field 


PM 


FIG. 2: The PPS encoding scheme for one input field is shown, where PNG denotes the pseudo¬ 
random number generator and PM denotes the phase modulator. 


modulated with a PPS AW The differenced signals of two output detectors are then summed 
over N and sampled to yield the decision variable. We can express the demodulation process 
in mathematical form 


/(A w ,A (r) ) 



k =1 


1, i — r 

0 , 


( 8 ) 


where Ajj' 1 , A^ are the PPSs of the input and reference fields respectively. The output 
decision variable is 1 if and only if A^, AjA'* are equal; otherwise the output decision variable 
is 0. The results are guaranteed by the properties of PPSs. If the input is a single-mode field, 
the scheme shown in Fig. [3] is employed to perform quadrature demodulation. Otherwise 
the scheme shown in Fig. [ 4 ] is used, in which the input field \ipi) = e* A<1> (cq |0) + /% |1)) is 
first splitted into two fields cqe* AW |0) and / 5je* A(!) 11), and two coherent detection processes 
are then performed on the two fields respectively. Noteworthily, the modes of the reference 
fields must be consistent with the two output fields. Thus there are two output decision 
variables a and /3, which correspond to the modes |0) and 11), respectively. We define ^a, 
as a mode status that displays not only existence of mode also phases due to the polarity 
of signals. Besides the quadrature demodulation, we can also easily measure the amplitudes 
of each modes cq|, |/3j| after mode spitting in the scheme. However, we will not deal with 
normalization of superposition coefficients a and (3 in our scheme, because the amplitudes 
of each modes are not important. 


II. SIMULATION OF MULTIPARTICLE QUANTUM STATES 


We discuss simulation of multiparticle quantum states using classical optical fields modu¬ 
lated with PPSs in this section. We first demonstrate that n classical optical fields modulated 






PNG A w 


a or p 



FIG. 3: The PPS quadrature demodulation scheme for one input field with single mode e* A |0) 
or e lX< ' l) |1) is shown. 



FIG. 4: The PPS quadrature demodulation scheme for one field with two orthogonal modes is 
shown, where the blue block MS denote the mode splitter. 

with n different PPSs can constitute a similar n2 n -dimensional Hilbert space that contains a 
tensor product structure [TO] . Besides, by performing quadrature demodulation scheme, we 
can obtain the mode status matrix of the simulating classical optical fields, based on which 
we propose a sequence permutation mechanism to simulate the quantum states. Classical 
simulations of some typical quantum states are then discussed, including product state, Bell 
states, GHZ state and W state. 

A. Classical optical fields modulated with pseudorandom phase sequences and 
their tensor product structure 

For convenience, here we consider two classical optical fields modulated with PPSs and 
their tensor product structure. Chosen any two PPSs of A^ and A® from the set S, any 
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two fields modulated with the PPSs can be expressed as follows, 

m = e iX(a) (cn a |0) + (5 a |1)), (9) 

\A) = e ix(b) (a b \0) + (3 b \l)). 

We define the inner product of two fields \ijj a ) and \ifj b ) as follows, 

N 

{^a\^b) = e ^ Afe _Afe ^ (“^a + ^o) ' ( 10 ) 

V k=1 

According to the properties of the PPSs, we can easily obtain 

f 1, a = b 

(^a\A) = { , ( 11 ) 

^ 0 , a 7 ^ b 

which shows that two fields modulated with two different PPSs are orthogonal. The or¬ 
thogonal property make contribution to the tensor product structure of multiple fields. The 
direct product states of the two fields can be expressed as follows, 

l^a) ® \A) = e< xW+xW ) ( a a a b |0> |0> + a a /3 b |0) |1) + Mb |1) |0) + Mb |1) |1)), (12) 

where A^ + A® remains in the set S due to the closure property. 

Quantum entanglement is only defined for Hilbert spaces that have a rigorous ten¬ 
sor product structure in terms of subsystems. As shown in Ref. ra. for any n classi¬ 
cal optical fields with n PPSs, the Basis for Hilbert space of simulation is spanned by 
|e* Ab) |iR 2 ... i n ) \j = 1 • • • **, ik = 0 or 1 j, with a total base state number of n2 n . It is ob¬ 
vious that the Hilbert simulation space is greater than what is required for simulation of 
quantum state. The tensor product structure and the efficient classical simulation of quan¬ 
tum entanglement have been discussed. Here we construct a similar structure of multiple 
classical optical fields based on the efficient classical simulation of quantum entanglement. 


B. Simulation of quantum states based on classical optical fields 


We have discussed quadrature demodulation process in Sec. |I C| Here we discuss how 
to simulate quantum state based on classical optical fields with the help of quadrature 
demodulation. 
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First, we consider the general form of n classical optical fields modulated with PPSs 
{A^P,..., A^} chosen from the set S, and the states can be expressed as follows, 


H’k) = ^ J a { k ) e lX{1 ' |o) + 




|i> 


(13) 


i=1 j=1 

It is noteworthy that although multiple PPSs are superimposed on both modes of the fields, 
all of the PPSs can be demodulated and discriminated by performing the quadrature de¬ 


modulation introduced in Sec. IC, which has already been verified by many actual commu¬ 
nication systems mm- 

Now we propose a scheme, as shown in Fig. [4j to perform the quadrature demodulation 


introduced in Sec. |I C[ In the scheme, quadrature demodulations are performed on each 
held, in which the reference PPSs are ergodic on |A^\ ..., A ^ }. Thus a mode status matrix 


M (^dq, /3? j , as shown in Fig. 5 can be obtained by performing n quadrature demodulations 
on the n classical optical fields. Of the matrix M (aj, each element is the mode status 
of the ith classical optical held when the reference PPS is A "A The element takes one of four 
possible values: (1,0), (0,1), (1,1) or 0, denote that the PPS A01 is modulated on mode 
|0) of the ith classical optical held, on mode |1), on both |0) and 11), on neither |0) nor 
|1), respectively. It is noteworthy that different modulation of the n classical optical helds 
correspond to different mode status matrixes, and vice versa. Thus we obtain a one-to- 
one correspondence relationship between the n classical optical helds and the mode status 
matrix. Besides, further discussion will show that structure of quantum states and quantum 
entanglement can be revealed in the mode status matrix, which means that a correspondence 
can also be obtained between the mode status matrix and quantum states. Thus we treat 
the mode status matrix as a bridge to connect the simulating helds and the quantum states. 

Now we consider how to construct the states based on M matrix, respectively. Thus we 
propose a sequence permutation mechanism to simulate each |T) based on M matrix, which 
is one of the simplest mechanisms for sequence ergodic ensemble. Assumed |T) contains n 
classical optical helds with n PPSs, namely the corresponding matrix contains n rows and 
n columns, the sequence permutation is arranged as 

Ri = {A (1) , A (2) ,..., A (n) } , R 2 = {A (2) , A (3) ,..., A (1) } ,..., R n = {A (n) , A (1) ,..., A (n - 1} } . 

(14) 

As shown in Fig. [6j we can obtain the mode status with same color for same sequence 
permutation, such as the red color corresponding to Ri, the blue color corresponding to i? 2 , 
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FIG. 5: The PPS quadrature demodulation scheme for multiple input fields is shown, where the 
DE block is shown in Fig. [4j 

etc. We obtain a direct product of n items for each R r , and the simulated quantum state is 
the superposition of the n product items. Therefore we obtain 

I*) = (a! |o> + ft 1 |i>) ®... ® (< |o) + k |i>) + (a? |o) + ft 2 |i>) ®® (ai |o) + ft |i)) 
+... + (a? |o> + ft |i>) ®... ® (a;- 1 |o> + ft- 1 |i>), (15) 

where IcftftJ is the mode status obtained from the matrix M ( a 2 . if’ 1 j. 

It is noteworthy that the mechanism we proposed above is one of the feasible ways to 
simulate the quantum state based on the mode status matrix. Other mechanisms may also 
work, as long as a sequence ergodic ensemble is obtained in the mechanism. The sequence 
permutation mechanism above can successfully simulate many quantum states, including the 
product states, Bell states, GHZ states and W states. We will discuss the related contents 
in next subsection. 


C. Simulations of several typical quantum states 

In this subsection, we discuss classical simulations of several typical quantum states, 
including product state, Bell states, GHZ state and W state. 
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A (n 1) 


A<» A® 


A (n) 



m 



+ 
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(«i|0> + ft 2 |l»®(a||0> + /? 2 3 |1» ...®(a*|0> + $|1» 

B'V 1 ■ 1 • , ’ i - 1 ■ 

+ b 

-|- “j- (a?|0> + ^ 1 n |l»®(ap0> + ^ 2 1 |l>)...®( 

aS -1 |0) + $J -1 |l>) 


FIG. 6: The scheme to simulate quantum state is shown, the mode status matrix M 
related to ith classical optical field and the reference PPS A Id, in which the mode status with same 
color for same sequence permutation. 

1. Product state 

First, we discuss classical simulation of n quantum product state. The simulation fields 
are shown as follows 

I’M = (|0> + |1» , (16) 


m = e“ w (|0) + |l». 


By employing the scheme as shown in Fig. [6j we obtain the mode status matrix 


/ 


M (aiJi) = 


( 1 , 1 ) 


\ 


( 1 , 1 ) ) 


(17) 
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which demonstrates that each classical optical field is the superposition of two orthogonal 
inodes and no entanglement is involved. According to Eq. (15), we obtain 


l*> = (|0> + |1»®...®(|0> + |1» (18) 

= | 0 ... 0 ) + | 0 ... 1 ) + ... + 11 ... 1 ), 

where \i\... i n ) = |A) 0 ... 0 | i n ), (i n — 0 or 1), which is same as a quantum product state 
expect a normalization factor. 


2. Bell states 


Now we discuss classical simulation of one of the four Bell states | v b + ) = 

(|0i> |O 2 ) + |li) 1 12 )), which contains two classical optical fields as follows 

|« = e ,A<1> |0)+e‘ xm |l), (19) 

life) = e u< 1 , | 0 }+e iJI 1 > |l). 


By employing the scheme as shown in Fig. [6j we obtain the mode status matrix 


m 


( 1 . 0 ) ( 0 , 1 )\ 

(0,1) (1,0) ) 


( 20 ) 


We note that in this case, the mode status matrix is irreducible, which corresponds to an 
entanglement state. According to the sequence permutation mechanism, we obtain that 
Ri = {A^b,A( 2 )} and R 2 = {A^A^}. Based on the mode status matrix, for the selection 
of Ri, we obtain |0) 0 10); for the selection of i? 2 , we obtain |1) 0 |1). If we randomly choose 
one selection between R\ and i? 2 , we can randomly obtain one result between |0) ® |0) and 
|1) ® 11), which is similar with the case of quantum measurement for the Bell state |'b + ). 
We can simulate the state except a norm based on the mode status matrix 


d' + ) — (aq |0) + f3\ |1) j <S> ^a 2 |0) + /3 2 |1) j + (a\ |0) + |1)^ 0 ^a 2 |0) + |1) j 


= |oo) + |n> 


( 21 ) 


which is same as the Bell state | v & + ) expect a normalization factor. 

In quantum mechanics, another Bell state | < h + ) can be obtained from |'1/ + ) by performing 
the unitary transformation a x : |0) -H- |1) on one of the particles. Using the same method, 
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we perform an unitary transformation on |^&) to flip its modes |0) |1). Thus we obtain 

two classical optical fields as follows 


iVh) = e a( 1 ) | 0 > + e a( 2 ) |l), ( 22 ) 

\^ 2 ) = e iA(2) | 1 ) + e a(1) | 0 ). 


By employing the scheme as shown in Fig. [ 6 j we obtain the mode status matrix 


M (a 1 ,,P () 


( 1 . 0 ) ( 0 , 1 ) \ 
( 1 . 0 ) ( 0 , 1 ) ) 


(23) 


According to the sequence permutation mechanism, here we obtain Ri and R 2 again. As 
the mode status matrix is different, for Ri, the result turns to be |0) g |1); for R 2 , we obtain 
|1) g |0). If we randomly choose one selection between R\ and R 2l we can also randomly 
obtain one result between | 0 ) g | 1 ) and | 0 ) g | 1 ), which is similar with the case of quantum 
measurement for the Bell state |<f> + ). We can simulate the state 


$ + ) — (a\ |0) + P\ |1) j g (aj |0) + P 2 |1)) + (“i |0) + Pi |1) j g (ot\ |0) + Pi |1)) 


= | 10 ) + | 01 ) 


(24) 


which is same as the Bell state |$ + ) expect a normalization factor. For other two Bell states 
IT - ) and they can be obtained from | v h + ) and |<f> + ) by using a 7 r phase transformation. 
We can distinguish |4/ _ ) and |$ _ ) from |T + ) and | < h + ) by using the signal polarity of 
quadrature demodulation. 


3. GHZ state 

For tripartite systems there are only two different classes of genuine tripartite entangle¬ 
ment, the GHZ class and the W class [ 21121 ]. First we discuss the classical simulation of 
GHZ state \^ghz) = ^ (|0i) |0 2 ) |0 3 ) + |li) 11 2 ) 1 13 )), which contains three classical optical 
fields as follows 

IV-.) = e““’|0)+e iJ 1 I ’|l>, (25) 

life) = e“‘”| 0 )+e iJ 13 > |l), 

life) = e iA< 1 ’| 0 )+e““ > |l). 
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Performing the scheme as shown in Fig. [6j we obtain the mode status matrix 




^(1,0) (0,1) 0 ^ 
0 (1,0) (0,1) 
v (0,l) 0 (1,0) , 


( 26 ) 


According to the sequence permutation mechanism, we obtain that Ri = (A 1 - 1 ), A^ 2 \ A®}, 
R 2 = {A^ 2 \ A^} and R 3 = {A^ 3 ), A^, A®}. Based on the mode status matrix, for the 

selection of Ri, we obtain |0) <£) |0) <g) |0); for the selection of R 2 , we obtain |1) <g) |1) <g) 11); for 
the selection of R%, we obtain nothing. Thus we can simulate the state based on the mode 
status matrix 


^GHz) - («} |0) + Pi |1)) <8> («2 |0) + I 1 )) ® («3 |0) + /3| |l)j (27) 

+ (&1 |0) + Pi |l)j ® (&2 |0) + @2 |1)) ® |0) + /3g |l)j 

= | 000 ) + | 111 ). 


4- W state 


Then we discuss the classical simulation of W state, 

I'M = A (111) |0 2 ) |0 3 ) + | 0l ) |1 2 ) |0 3 > + |00 |0 2 > |1 3 » , (28) 

which contains three classical optical fields as follows 

iVh) = e a(1) |1) + e ix(2) |0) + e ix(3) |0), (29) 

|^ 2 > = e a(1) |1) + e ix(2) | 0 ) + e ix(3) | 0 ), 

|^ 3 > = e a(1) |1) + e ix(2) | 0 ) + e ix(3) | 0 ). 


It is noteworthy that the three classical optical fields can be produced from one single field by 
using two beam splitters, which is quite similar with the generation of W state in quantum 
mechanics. Performing the same scheme, we obtain the mode status matrix 


M(y,/3f) 


' (o,1) (1,0) (1,0)^ 
(0,1) (1,0) (1,0) 
v (0,l) (1,0) (1,0) j 
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(30) 






According to the sequence permutation mechanism, we use R \, R 2 and R 3 again. Based on 
the mode status matrix, we obtain |1)<8> |0)<E> |0), |0) <8> |0) (8) 11), |0) 011) 0 |0) for the selection 
of i?i, i? 2 , i? 3 , respectively. We find an interesting fact that when we obtain the state |1) 
of the first field, R\ must be selected, thus only the |0) ® |0) state can be obtained from 
the other two fields; otherwise when we obtain the state |0) of the first field, the selection 
can be R 2 or R 3 , thus the state of |0) <8> |1) + |1) <8> |0) can be obtained from the other two 
fields. This fact is quite similar with the case of quantum measurement and the collapse 
phenomenon for W state in quantum mechanics. We can simulate the state based on the 
mode status matrix expect a normalization factor, 

l^w) = |0> + Pi |1)) O (ctl |0) + Pl |1)) ® |0) + Pi |1)) (31) 

+ (“1 |0) + Pi |1)) ® |0) + Pi |1)) ® (“3 |0) + Pi |l)j 

+ (“1 |0> + Pi |l)j ® {ol\ |0) + Pi |l)j ® (®3 |0) + Pi |l)j 
= | 100 )+ | 010 ) + | 001 ). 

Above we discussed the possibility of simulating several typical quantum states by clas¬ 
sical fields and mechanism of coherent detection. 

III. SIMULATION OF QUANTUM COMPUTATION 

In this chapter we will propose a method to simulate quantum computation. In quantum 
computation, any quantum state can be obtained from initial states by using unitary trans¬ 
formation of universal gates. Similarly, we can construct simulation of all kinds of quantum 
states by using a gate array, such as GHZ state and W state, even to realize Shor’s algorithm 
and Grover’s algorithm. 

A. Gate array model to simulate quantum computation 

In Ref. [TO], a constructure pathway of simulation states is shown. Here we use the 
same model as shown Fig. [Tj however a gate array (GA) is proposed instead of the unitary 
transformation to produce simulation states. Now we discuss some basic units of gate array 
model. 

(1) Combiner and splitter 


17 


’e' 2 ® \0)+Ifi? e liW |1) 


|» t >=e M(1) |0> 

I02>=e“*l<» 


l0n>=e i ' i(,,) |O> 


e M(1) (|0>+|l)) 

e u(2) (|0>+|D) 


e “ (n) (|0>+|D) 



l^2)=Z«2 )e “ ( ° |0)+^ ) e u( ' ) |1> 


|0>+Z/»?e uW ID 


FIG. 7: A gate array model to simulate quantum computation is shown, where GA denotes the 
gate array. 

Different from quantum state, we can conveniently combine and split a classical field 
by using an optical splitter device. Therefore we define two basic device as combiner and 
splitter as shown in Fig. [8] (a) and (b), respectively. 


w 



(a) 

l*<> 

M) = i<m 

—«-» 


W) = iV’i) 


(b) 


FIG. 8: Two basic devices (a) combiner and (b) splitter are shown. 


(b) Mode control Gates 

Further, we define 4 kinds of mode control gates as selective mode transit devices with 
one input and one output as shown in Fig. [9] They are defined by following: 


GateA : |V’i) — e lX (|0) + |1)) —> |V0 — 0, (32) 

GateB : \A) = e iXl " (|0) + |1» -> |</>') = e ttW |0>, 

GateC : |*) = e“ M (|0) + |1» -> |$) = e“ w |1), 

GateB : |*) = e““’ (|0) + |1» -> |$) = e“ 0> (|0) + |1». 


Now we will discuss a basic structure of gate array model. According to Sec. IIC 2 
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Gate A 


? uC0 (|0)+|l» 


WM 


■ 


! uC0 (|0>+|l)) 

a 

W)=e u<0 |0> 

! “ C0 (|0>+|1» 

■ 

W)= e ‘ 1<0 |l> 


Hf 


A°(io)+ii» 

—H- 

WH uCO (lo)+|i» 


Gate B 



FIG. 9: Mode control gates as selective mode transit devices are shown. 


we can simulation of quantum states by using sequence permutations. Similar to field 
programmable gate array (FPGA), we propose a simple structure of gate array to obtain 
permutation structure as shown in Fig. 10 Gates Gkj constituted by the basic units can 


transform 1^) to achieve certain |?/4). It is easy to know that a sequence permutation with 
circulation of p + 1 needs at least p combiner devices and 2 p control gates. Any states is 
capable to be constructed by applying this structure that will be strictly proved in future 
paper. 


I'M 


Wk+1) mod n) 





l$(fc+2) mod n) 



* W) 


I^Cfc+p) tnodn) 


FIG. 10: A basic structure of gate array model is shown. 


Finally, we illustrate two gate array models to transform product states to GHZ state 


and W state as shown in Fig. 11 


B. Simulation of quantum algorithm 

1. Shor's Algorithm 

Shor’s algorithm is a quantum algorithm for integer factorization that runs only in poly¬ 
nomial time on a quantum computer [T%1 IT*?] . Specifically it takes time and quantum gates 


19 


























(A) 



(B) 


FIG. 11: Gate array models to transform product states to (a) GHZ state and (b) W state are 
shown. 

of order 0((logiV) 2 (loglogiV)(logloglogiV)) using fast multiplication, demonstrating that 
the integer factorization problem can be efficiently solved on a quantum computer and is 
thus in the complexity class BQP. 

In this chapter, we discuss the simulation of Shor’s algorithm. A algorithm similar Shor’s 
algorithm is demonstrated, that factored N = 15 into 3x5, using the gate array model with 
8 classical optical fields. First, we chose a random number a coprime with 15, for example 
a = 7. We define a function as followed 


/ (x) = a x modN = 7 x modl5. 


(33) 


The key step of the Shor’s algorithm is to obtain the period r to satisfy 


f (x + r ) = 7 x+r modlh = 7 x modlh — f (x). 


(34) 


In order to construct / (x), we prepare 8 classical optical fields modulated with 8 PPSs: 


l«=e ,A “’(|0) + |l»,fc=1...8. 


(35) 


We can express the product state as followed: 


255 


I'p) = hM®|^ 8 > = E tf> • 

3=0 


(36) 


Further, we construct a gate array model as shown in Fig. 12 After passing through the 


gate array and the classical fields will become the forms as follows: 

K) = fe a(1) +e' A(2) +e a(3) +e a(4) )(|0) + |l)), 


(37) 
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W 2 ) = (e a(2) +e iA(3) +e iA(4) +e iA(5) )(|0) + |l)), 

|^ 3 ) = ( eU(3) + eiA(4) ) 1°) + ( e ' A(5) + eiA(6) ) I 1 ), 

m = (e^ C4) + e iA(6) ) |0) + (e^ A(5) + e iA(7) ) |1), 
|^) = (e a(5) + e ix(6) + e iA(7) ) |0) + e iA(8) |1), 
l^e) = e U(6) |0) + (e iA(7) + e a(8) + e* A(1) ) |1), 

\if/ 7 ) = (e iA(7) + e iA(1) + e iA(2) ) |0) + e iA(8) |1), 

\i/Q = (e ut8) + e iA(2) ) |0) + (e a(1) + e‘ A(3) ) |1). 


By using the quadrature demodulation, we can obtain the mode status matrix: 




( ( 1 , 1 ) ( 1 , 1 ) ( 1 , 1 ) ( 1 , 1 ) 0 0 0 

0 ( 1 , 1 ) ( 1 , 1 ) ( 1 , 1 ) ( 1 , 1 ) 0 0 

0 0 ( 1 , 0 ) ( 1 , 0 ) ( 0 , 1 ) ( 0 , 1 ) 0 


0 0 0 

0 0 0 

(0,1) 0 0 

(1,0) (1,0) 0 

V (0,1) (1,0) (0,1) 


( 1 , 0 ) ( 0 , 1 ) ( 1 , 0 ) ( 0 , 1 ) 


0 

0 

0 

0 


\ 


0 ( 1 , 0 ) ( 1 , 0 ) ( 1 , 0 ) ( 0 , 1 ) 

0 0 ( 1 , 0 ) ( 0 , 1 ) ( 0 , 1 ) 

0 0 0 ( 1 , 0 ) ( 0 , 1 ) 

0 0 0 0 (1,0) y 


(38) 


Using the scheme mentioned in Sec. IIB we can obtain the simulated states: 


*') = (|0) + 14) + 18) + 112)) |1) (39) 

+ (|1) + |5) + |9) + 113)) |7) 

+ (|2) + |6) + |10) + |14))|4) 

+ (|3) + 17) + 111) + 115)) 113). 


There are four kinds of superposition classified from last four qubits containing the values 
of/ (x) (|1), |7), |4) and 113)) in output states, which means the period of / (x) = 7 x modl5 
is r = 4. It is worth noting that, different from quantum computing, we obtain the expected 
period of without operating quantum Fourier transformation. The remaining task is much 
easier. Because N = 15, a = 7, r = 4, we obtain 

gcd (a 5 ± 1, N) = gcd (7 5 ± 1,15) = gcd (49 ± 1,15), (40) 
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where gcd(48,15) = 3, and gcd(50,15) = 5. Finally, we can deduced that N (15) = 3x5. 

Here, this is the simulation of quantum Shor’s algorithm. Apparently, this algorithm 
obtains the period of / (x) without quantum fourier transform and takes time and gates 
of order 0((logAt) 2 ) by rough estimation. We will discuss this problem in detail in future 
paper. 


|t/. 1 >=e i ' ,tl) (|0>+|l» 

l</- 2 >=e i,,(2) (io>+|i» 

l<M=e i ' ,(4) (|0>+|l» 

|^ 5 )=e i ' l(S) (|0>+|l» 

|^ 6 >=e i - ,<6) (|0>+|l)) 

IV'7>=e a(7) (|0>+|l» 

IV. s >=e i ' ,t8) (|0>+|l» 



|i/>i)=(e“ (1) +e“ t2) +e uC3) +e uC4, )(|0>+|l» 
| xj,’ 2 )=(e^ m +e iiW +e aW +c“ (5) ) (| 0)+11» 
|,/,')=(e i;|(3) +e“ (4) )|0>+(e i ' ,tS) +e i;it6) )|l) 
|^;)=(e u(4) +e“ W )|0)+(e l;l(5) +e u(7) )|l) 

|^)=( e “ (5 He“ C6 Ae u(7) )|0>+ e “ (8) |l> 

|^)=e u(6) |0>+( e “ (7) +e “ C8) +e u(1) )|l> 

|Vi^=(e u(7) +e uti:l +e uC2) )|0>+e u(8) |l) 

+e uC2) ) |0)+(e u(1) +e “ (3, )n> 


FIG. 12: A gate array model to realize Shor’s algorithm is shown. 


2. Grover’s Algorithm 

Grover’s algorithm is a quantum algorithm for searching an unsorted database with N en¬ 
tries in 0(N 1 ^ 2 ) time and using 0(log N) storage space [20]. Grover’s algorithm demostrates 
that in the quantum model searching can be done faster than classical computation; in fact 
its time complexity 0(N 1 ^ 2 ) is asymptotically the fastest possible for searching an unsorted 
database in the linear quantum model. However, it only provides a quadratic speedup rather 
than exponential speedup over their classical counterparts. 

We now discuss classical simulation of Grover’s algorithm. According to Grover’s algo¬ 
rithm, the steps of the simulation are given as follows. Let | S) denote the uniform superpo¬ 
sition over N states, 

N 

is) = y , w < 41 ) 

i —1 
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where \xi) are the states. For example, we obtain \S) as a superposition state of 13 random 
numbers: 


\S) = |61) + 163) + 1117) + 1125) + 1140) + 1142) + |148) (42) 

+ 1212) + 1187) + 159) + 1238) + |247) + |76) 

= 100111101 )+ 100111111 ) + 101110101 )+ 101111101 ) 

+ 110001100 ) + 110001110 ) + 110010100 ) + 111010100 ) 

+ 110111011 ) + 100111011 ) + 111101110 ) + 111110111 ) 

+ 101001100 ). 


We choose 8 classical fields modulated with 8 PPSs and after passing through a suitable 
gate array, that become the forms as follows: 


K> 

m 

m 

m 

m 

m 

IV4> 


e ixW + e u(2) + + A(5) + e^ (8) l |0) + (+ A(3) + + A(4) + + A(5) + + A(6) + e^') |1), (43) 


= e 




+ e u(4) + e a(5) + e 


a(°) 


|0) + (e iA(1) + e lAW + e* AW + e lA '" + e lAW ) |1), 


,a( 3 ) 


A(S) , 




= e 




+ e‘ A(5) + e a(6) 


|0) + (e aU) + e* A(3) + e iAW + + e 4A '°' ) |1), 


i\W 


i\W . J\(S) 


= e 


+ e + e* 


+ e l 


|0)+(e iA(2) +e lAW +e lAW +e jA "M|l), 


i\W 


A(5) 


i\W 


e ix(3) + e iA(6) + e iA(8) ) |0) + (e* A(1) + e lA(2) + e* A(4) + e‘ A(5) + e‘ A(6) + e * A(7) ) |1), 
e a(2) |0) + (V A(1) + e a(3) + e iA(4) + e a(5) + e a(6) + e a(7) + e a(8) ) |1), 
e iA(1) + e iA(2) + e u(6) + + A(?) + e^) |0) + (e ixW + e a(3) + + A(4) + + A(5) + + A(7) 


|1), 


= e 


+ e + e l 


+ e l 


|0)+ (e* A + e* 


i\(e) , iA(S) 

+ e + e 


| 1 )- 


Due to the scheme mentioned in Sec. |II B[ any simulated state must correspond to a 
certain sequence permutation. Therefore, the problem to determine whether |x) exists in 
\S) become that to search the corresponding sequence permutation. For example, we search 
the number |x) = 1148) = 110010100) in (S'). First, the classical fields of state \S) pass 
through a gate array as shown in Fig. [13] Then we obtain the mode status matrix by using 
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the quadrature demodulation: 


Mfap}) 


0 

0 

( 0 , 1 ) 

( 0 . 1 ) 

(0 

0 

( 1 . 0 ) 

0 

( 1 , 0 ) 

(1 

0 

( 1 . 0 ) 

0 

0 

(1 

0 

( 0 . 1 ) 

0 

( 0 . 1 ) 

(0 

0 

0 

( 1 . 0 ) 

0 


>, 1 ) 

0 

( 0 . 1 ) 

( 0 . 1 ) 

(0 

, o ) 

( 1 . 0 ) 

0 

0 


0 

( 1 . 0 ) 

( 1 . 0 ) 

0 

(1 


, 1 ) ( 0 , 1 ) ( 0 , 1 ) 0 \ 
, 0 ) ( 1 , 0 ) 0 0 
, 0 ) ( 1 , 0 ) 0 0 
, 1 ) 0 ( 0 , 1 ) 0 
0 ( 1 , 0 ) 0 ( 1 , 0 ) 
, 1 ) ( 0 , 1 ) ( 0 , 1 ) ( 0 , 1 ) 

0 ( 1 , 0 ) ( 1 , 0 ) ( 1 , 0 ) 
, 0 ) 0 ( 1 , 0 ) 0 j 


(44) 


Finally, it is easy to search only corresponding sequence permutation R 4 = 



If we 

choose x) 

= 240) = |11110000) : 

obtain the mode status matrix 








( 0 

0 

(0.1) 

(0,1) 

(0.1) 

( 0 , 1 ) 

(0,1) 0 > 



(0,1) 

0 

(0.1) 

0 

(0.1) 

0 

(0,1) (0,1) 



(0,1) 

0 

(0.1) 

(0,1) 

0 

0 

(0,1) (0,1) 


M(aiJi) = 

0 

(0,1) 

0 

(0,1) 

(0.1) 

0 

(0,1) 0 



0 

0 

(1,0) 

0 

0 

(1,0) 

0 (1,0) 



0 

(1,0) 

0 

0 

0 

0 

0 0 



(1,0) 

(1,0) 

0 

0 

0 

(1,0) 

(1.0) (1,0) 



l 0 

(1.0) (1.0) 

0 

(1,0) 

0 

(1,0) 0 ) 



we can 


(45) 


In the mode status matrix, we can not search any corresponding sequence permutation. 
Therefore we can conclude |x) = |240) does not exist in I5 1 ). 

Here, this is the simulation of quantum Grover’s algorithm. Different from Grover’s 
algorithm, this algorithm for searching an unsorted database with N entries in 0((loglV) 2 ) 
time and using 0((log N) 2 ) storage space by rough estimation. We will discuss this problem 
in detail also in future paper. 


IV. CONCLUSIONS 

In this paper, we have discussed a new scheme to simulate quantum states by using classi¬ 
cal optical holds modulated with pseudorandom phase sequences. We first demonstrated that 
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fe' iC2) +e u(5) +e u(6 ) 
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r e a(» +e u(«) +e u(5) +e uP) 
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f e u 1 2 (3) + e“ (6, +e“ <8) 
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( e uW +e u«> +e u») +e uW +e u«) 


)| 0 > 


(e u<2) +c u<3) +e u<5) +e u<7) )|0) 




FIG. 13: A gate array model to select \x) from | S) is shown. 

n classical optical fields modulated with n different PPSs can constitute a n2”-dimensional 
Hilbert space that contains tensor product structure similar to quantum systems. Further, 
by performing quadrature demodulation scheme, we obtained the mode status matrix of 
the simulating classical optical fields, based on which we proposed a sequence permutation 
mechanism to simulate the simulated quantum states. Besides, classical simulation of several 
typical quantum states was discussed, including product state, Bell states, GHZ state and W 
state. Finally, we generalized our simulation and discussed a generalized gate array model to 
simulate quantum computation. The research on simulation of quantum states is important, 
for it not only provides useful insights into fundamental features of quantum mechanics, but 
also yields new insights into quantum computation and quantum communication. 
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Fig. 1 The time sequence relationship of the PPS is shown. 


Fig. 2 The PPS encoding scheme for one input field is shown, where PNG denotes the 
pseudorandom number generator and PM denotes the phase modulator. 

Fig. 3 The PPS quadrature demodulation scheme for one input held with single mode 

e lX<l) |0) or e* A(l |1) is shown. 

Fig. 4 The PPS quadrature demodulation scheme for one held with two orthogonal modes 
is shown, where the blue block MS denote the mode splitter. 

Fig. 5 The PPS quadrature demodulation scheme for multiple input helds is shown, where 
the DE block is shown in Fig. |4} 

Fig. 6 The scheme to simulate quantum state is shown, the mode status matrix M /d^j 
related to zth classical optical held and the reference PPS in which the mode status 
with same color for same sequence permutation. 

Fig. 7 A gate array model to simulate quantum computation is shown, where GA denotes 
the gate array. 

Fig. 8 Two basic devices (a) combiner and (b) splitter are shown. 

Fig. 9 Mode control gates as selective mode transit devices are shown. 

Fig. 10 A basic structure of gate array model is shown. 

Fig. 11 Gate array models to transform product states to (a) GHZ state and (b) W state 
are shown. 

Fig. 12 A gate array model to realize Shor’s algorithm is shown. 

Fig. 13 A gate array model to select |x) from \S) is shown. 


